Two-Dimensional Melting of a Crystal of Ferrofluid Spikes 
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We report the observation of the transition from an ordered solid-like phase to a disordered liquid- 
like phase of a lattice of spikes on a ferrofluid surface submitted to horizontal sinusoidal vibrations. 
The melting transition occurs for a critical spike displacement which is experimentally found to follow 
the Lindemann criterion, for two different lattice topologies (hexagonal and square) and over a wide 
range of lattice wavelengths. An intermediate hexatic-like phase between the solid and isotropic 
liquid phases is also observed and characterized by standard correlation functions. This dissipative 
out-of-equilibrium system exhibits strong similarities with 2D melting in solid-state physics. 

PACS numbers: 47.65.Cb, 64.70.D-, 05.70.Ln 



Melting of 3D crystals has been a substantial field of 
interest in condensed matter physics for over a century 
[l|. In 1910, Lindemann assumed that melting occurs 
when the thermal vibrations of atoms make them collide 
with each other [2.]. Then, it was suggested that the crit- 
ical rms value of the atomic displacement at the melting 
transition is rather given by a constant fraction {c^ 10 %) 
of the interatomic distance Q. This assumption, now 
known as the Lindemann criterion, has been succesfully 
used to predict melting temperatures of a wide range of 
3D crystals [3] . On the other hand, melting of a 2D solid 
is a less understood phenomenon since 2D lattices do not 
display true long-range translational order at finite tem- 
peratures [s*] . Consequently, the Lindemann criterion was 
then long thought inapplicable in 2D [6]. A 2D melting 
theory driven by topological defects has then been de- 
veloped by Kosterlitz, Thouless, Halperin, Nelson, and 
Young (KTHNY) 0] • Numerical studies have then shown 
the validity of the Lindemann criterion in 2D [6] , and con- 
sistency with the KTHNY scenario [8]. Although most 
experimental observations are consistent with this the- 
ory, it is quite difficult to establish unambiguously the 
existence of second-order solid-to-hexatic-to-liquid tran- 
sitions [9]. Consequently, the melting transition in 2D 
systems remains an active research field in solid-state 
physics [10], and in numerous domains including thin 
colloidal suspensions 
granular monolayers 



12|, liquid films [13], vibrated 
15 1, magnetic solid films [l6|. 



and vortex lattices in superconductors [l7| . 

A ferrofluid is a stable suspension of nanometric mag- 
netic particles diluted in a carrier liquid which displays 
striking phenomena such as flows driven by a mag- 
netic field gradient, magnetic levitation, labyrinthine and 
Rosensweig instabilities [18]. This latter occurs when a 
normal static magnetic field applied to a pool of ferrofluid 
exceeds a critical value: the flat free surface becomes 
unstable and a stationary hexagonal pattern of surface 
spikes grows. Following the pioneer work of Bragg et al. 
for the assemblage of soap bubbles [l9] , one can consider 



the Rosensweig spike lattice as a macroscopic analogous 
of a 2D crystalline structure and expect solid-like be- 
haviors despite the complexity of the ferrohydrodynamic 
interaction between spikes and the dissipative nature of 
the lattice. An interesting feature of this system is that 
both the lattice wavelength and topology can be tuned 
by a single external control parameter (see below). 

In this Letter, we report the first observation of a tran- 
sition (melting) between an ordered phase (solid) and a 
disordered phase (liquid) of a lattice of spikes on the sur- 
face of a ferrofluid submitted to sinusoidal vibrations. 
We study which parameter controls the transition and 
whether a Lindemann criterion can be applied. We char- 
acterize structural changes across this transition using 
classical condensed matter physics concepts. Our system 
being dissipative with nonequilibrium steady states, the 
comparison of 2D melting in out-of-equilibrium systems 
and in equilibrium ones is of primary interest [13, llH . 

The experimental setup has been described previously 
0. It consists of a container filled with a ferrofluid up 
to a depth h = 2 cm. In order to discriminate any finite 
size or boundary condition effects, containers of differ- 
ent shapes and sizes are used: cylindrical containers i) 
20 cm or ii) 12 cm in inner diameter, and iii) a rectan- 
gular container 13 x 9 cm^ sides. All containers are 4 
cm depth. The ferrofluid used is a aqueous suspension 
of maghemite particles [20]. Its properties are: density, 
p = 1324 kg/m^, surface tension, 7 = 59 x 10~^ N/m, 
initial magnetic susceptibility, Xi = 0.69, magnetic satu- 
ration Msat = 16.9 X 10^ A/m, the viscosity being close 
to the water one. The container is placed in a vertical 
magnetic induction B generated by two horizontal coax- 
ial coils (up to 780G) [2u\. A pattern of spikes on the 
surface of the ferrofluid is observed when B is above a 
critical value Be = 294 ± 2 G. This is close to the the- 
oretical value of 292.3 G computed as the threshold of 
the Rosensweig instability of our ferrofluid 0]. We de- 
note the dimensionless magnetic induction 5* = B/Bc. 
Values 5* > 1 will be used in the following. A hexagon- 
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square transition occurs at a second threshold [2l|, |22| , 
and is observed here at = 1.45. For a fixed 5*, the 
lattice of spikes of ferrofiuid is vibrated by means of the 
horizontal motion of a rectangular Tefion plate (plung- 
ing perpendicularly to the fiuid at rest) that is driven 
sinusoidally by an electromagnetic vibration exciter at 
frequency / and maximal displacement amplitude a in 
the ranges 5 < / < 50 Hz and 0.5 < a < 7 mm, re- 
spectively. The vibrating plate acceleration is measured 
with an accelerometer. A high resolution camera located 
above the ferrofiuid container allows us to measure the 
spike positions by calculating centers of the bright spots 
produced by refiections of light on the top of each spike 
(see Fig. [T] for a typical snapshot). A particle tracking 
method is used to measure the spike displacement d with 
respect to the vibrating plate displacement a. For the 
three containers used, d is found to be a linear function 
of a with a frequency-independent coefficient a = d/a: i) 
a = 0.77 ± 0.04; ii) a = 0.85 ± 0.05; iii) a = 0.65 ± 0.05. 
The slight variations of a with the container may be due 
to different relative positions of the vibrating plate to the 
front row of spikes. Since the vibrations are sinusoidal, 
the acceleration amplitude of the vibrating plate then 
reads F = 47r/^a = Aupd/a. 

Figure [T] shows photos of the spike lattice on the fer- 
roffuid surface for two different vibrating plate accelera- 
tions F, at a fixed B* . At smah F (Fig. [1^), the lattice 
vibrates at the same frequency / than the vibrating plate 
one: each spike displacement fiuctuates about its equilib- 
rium position as would do atoms of a crystal due to ther- 
mal motion. When F is enough increased (at a constant 
/), the ordered lattice melts: ferrofiuid spikes do not have 
a steady position anymore and rows of spikes slide past 
each other. At still higher F (Fig. [TJd), the spikes display 
highly disordered dynamics and two (or more) colliding 
spikes can fuse into one larger unstable spike. 

The onset of melting is defined once a spike has a mo- 
tion larger than the lattice wavelength. It appears to be a 
sharp transition: the corresponding critical acceleration 
at the melting F^ can be measured with an experimental 
error of less than 5%. No hysteretic behavior is observed. 
Figure [2] then shows the evolution of the critical acceler- 
ation F^ as a function of the vibration frequency / for 
5 different applied magnetic inductions. For 8 < / < 20 
Hz, Tm is found to follow a power law with a scaling 
exponent 2 ± 0.1. This result means that along the melt- 
ing transition, the quantity F^//^ is a constant. Since 
the forcing is sinusoidal, this means that the amplitude 
of the spike displacement at melting, d^, is the relevant 
parameter for the transition: 
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For each 5*, the value of dm is extracted, using Eq. ([T]), 
from the ordinate intercept of the slope of each log-log 
plot in Fig. [21 The inset of Fig. [2] then shows the depen- 



FIG. 1: Top views of the spike lattice on the surface of a 
ferrofiuid for two different sinusoidal forcing amplitude: (a) 
hexagonal solid-like phase (F = 3 m.s~^), 5 ^ 15.8 mm, (b) 
liquid-like phase (F = 20 m.s~^). The vibrating plate is shown 
in white, on the left-hand side. / = 8 Hz, B* — 1.2. 



dence of the critical displacement dm (rms value) with the 
dimensionless magnetic induction B* for different con- 
tainers, dm is found to increase with 5* from 1.6 mm 
to 2.6 mm. Note that these values are one order of mag- 
nitude lower than the typical wavelength between spikes 
Ac = 27rY^7/ {pq) 13.4 mm at the Rosensweig instabil- 
ity threshold [18]. 

In order to establish a criterion for the melting transi- 
tion, both the lattice wavelength. A, and the spike height, 
/i, are measured to be compared to the spike displacement 
at melting, dm- The amplitude h of the ferroffuid spike 
is measured by means of a capacitive wire gauge |20| . 
Right inset of Fig. [3] shows as a function of 5*: data 
can be well fitted by — 1 (see dashed line) in good 
agreement with theory [21] and a previous observation 
[2^. The melting displacement, dmiB"") (inset of Fig. 
[2]), clearly does not follow the same law than the spike 
amplitude one, h{B^) (right inset of Fig. [3|), both quali- 
tatively and quantitatively. This means that the melting 
mechanism is not induced by surface waves higher than 
spike height. For a fixed 5*, the wavelength A is deter- 
mined by measuring the average spike-to-spike distance 
6 on experimental frames (see Fig. [1]). From simple ge- 
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FIG. 2: (color online). Critical acceleration Fm at the melting 
as a function of the vibration frequency / for different dimen- 
sionless magnetic induction B*: (x) 1.09, (+) 1.39, (o) 1.52, 
(□) 1.77 and (o) 2.03. Dashed lines have slopes 2 (constant 
displacement). Inset: Critical displacement dm as a function 
of B* for 3 different containers: (o) i, (x) ii, and (□) iii. 
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FIG. 3: (color online). Lindemann ratio as a function 
of the dimensionless magnetic induction B* measured for 3 
different containers: (o) i, (x) ii, and (□) iii. The dotted 
horizontal line is 7m = 0.14. Left inset: Lattice wavelength, 
A, as a function of B*. Right inset: spike amplitude, h, as a 
function of B* . \/B* — 1 fit (dashed line). 



ometrical considerations, one has A = {\^/2)S for the 
hexagonal pattern {1 < B* < 1.45), and X = 6 for the 
square pattern (5* > 1.45). Left inset of Fig. [3] then 
shows A as a function of B*. A is found to increase with 
B* as previously reported [22[. Let us now define an 
analog of the Lindemann ratio in solid-state physics. 



X{B*) 



(2) 



that is the ratio between the rms value of the spike dis- 
placement, dm^ at the melting transition and the lattice 
wavelength A at a fixed 5*. As shown in Fig. [3l 7^ 
is found to be independent of 5*, even at the hexagon- 
square transition. This is the first experimental obser- 
vation of the Lindemann criterion for the 2D melting 
transition of a crystal of ferrofluid spikes. The melting 
Lindemann ratio 7^ is found to be equal to 0.14 ± 0.02 
which is in the range of 3D crystalline solid values (0.1 
- 0.2) [1] and close to the reported value in 2D granular 
fluids (0.15) Note that defining the Lindemann ra- 
tio with respect to the lattice wavelength. A, instead of 
the spike-to-spike distance, leads to a single value for 
both lattice topologies (hexagonal and square). 

To characterize structural changes along the melting 
transition at a fixed 5*, one computes standard po- 
sitional and orientational correlation functions, respec- 
tively, g{r) = {n{r')n{r + r'))/(n(r'))^, where n is the 
spike density at a distance r from a reference (brack- 
ets are an average over the spatial variable r'), and 
gQ{r) = (^g(r')^6(^ + where the star denotes 

the complex conjugate and ^6(^j) = {ex.p{i66jk))k with 
6jk the orientation angle of the bond between the centers 
of the spike j and of the neighbouring spike k [9|] . These 
functions are averaged over 200 frames leading to an error 



bar of 2%. Figure [4] shows experimental curves of the ra- 
dial distribution functions g{r) for different values of the 
dimensionless forcing parameter e = — jrn) / Im where 
J = d/X. For e < (before the melting), g{r/S) displays 
characteristic features of a hexagonal structure: the first 
spikes positions are in very good agreement with values 
predicted from simple geometrical calculations r/S = 1, 
>/3, 2, >/7, 3 and >/T2 (see dashed lines in Fig. [4^). For 
£ = 0, the positional order is clearly short-ranged al- 
though some of the characteristic spikes of the hexagonal 
lattice remain visible: r/S = 1, V^, y/7 and a/T2 (see 
dashed lines in Fig. [4)3). As e is further increased, the 
positional-order range becomes shorter, and for 6: = 1, 
only the characteristic spikes of an isotropic liquid phase 
are observed: r/6 = 1, 2 and 3 (see dashed lines in Fig. 
Ht). These typical structural changes show strong sim- 
ilarities with the ones reported during numerical simu- 
lation of hard-disk fluid when the 2D solid-liquid phase 
transition is approached [13]. Our results are also con- 
sistent with the KTHNY theory in solid physics which 
predicts the existence of a hexatic phase between the 
crystalline and the liquid phases, characterized by a long- 
range orientational order (algebraic decay) and a short- 
range positional order (exponential decay) [7]. Indeed, 
right-hand side insets of Fig. [H show the orientational 
correlation function, g6{r). Just above the melting tran- 
sition {e > 0), the long-range orientational order (well 
fitted by an algebraic decay) is preserved as required for 
a hexatic phase (see Fig. HJd). Note that to utterly dis- 
criminate it from an exponential decay, a larger number 
of lattice periods should be necessary [IJ]. For £ = 0, 
one has geir) ~ 7--0-20 agreement with KTHNY sce- 
nario as the power-law exponent is predicted to continu- 
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FIG. 4: (color online). Radial distribution functions g{r/5) 
for different forcing: (a) e — —0.4, (b) s = 0, (c) £ = 1. 
B* — 1.1. f — 10 Hz. 6 — 15.5 mm. Dashed lines show 
predicted maximum positions of g{r) (see text). Right insets: 
Log- log orientational correlation functions gQ{r/S). Dashed 
lines are (a) 0.89, (b) (c) exp[-0.65(r/^)]. Left 

insets: typical thresholded photos of the ferrofluid surface. 



ously decrease in the hexatic phase down to —0.25 at the 
hexatic-liquid transition As e is further increased, the 
decay rate increases strongly, and for £ = 1, the orienta- 
tional order is short-ranged (well fitted by an exponential 
decay) as expected for a liquid phase. 

Such an intermediate hexatic phase between the crys- 
tal and liquid ones has been also reported in 2D melting 
of colloidal crystals with absolutely calibrated interaction 
in good agreement with the KTHNY theory [IJj]. Note 
that the 2D melting of these systems can depend on the 
interaction with their carrier substrate Our results 
put forward that both the Lindemann criterion and the 
KTHNY scenario are applicable to a more complex 2D 
system such as a crystal of ferrofluid spikes with com- 
plex ferrohydrodynamic interaction. Such a continuous 
solid-liquid transition, via a hexatic phase, also strongly 
differs of the first-order transition in 3D systems. Finally, 
our work emphasizes the analogy between 2D melting of 
equilibrium systems and out-of-equilibrium steady state 
ones. Such a correspondence has been observed in a 2D 
granular fluid IJ, [l5| , where spatial homogeneity of the 
energy injection within the system is underlined to be 



the main ingredient for these similarities with equilibrium 
dynamics [l5i]. Our results suggest that equilibrium-like 
properties can be observed even though energy injection 
is inhomogeneous. This should deserve more studies to 
have a complete description of the 2D melting transition 
of a dissipative crystal. 
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